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Classical transport theory for colored particles is investigated and employed 
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' our approximation scheme. We then derive the generating functional of hard 

c~| ' thermal loops from a constraint satisfied at leading-order by the color current. 
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I. INTRODUCTION 



Following a recent Letter [jT|, we present an expanded and self-contained account of 
the derivation of the hard thermal loops (HTLs) of QCD from classical transport theory. In 
addition, we justify the use of the ad hoc phase-space integration measure for classical colored 
particles. This justification is based on the phase-space symplectic structure, and relates 
directly the (dependent) color charges to a set of (independent) Darboux variables. We also 
discuss formally the gauge invariance properties of the system of coupled non-Abelian Vlasov 
equations, and exploit the gauge principle to justify the approximation scheme we use. In 
order to show how physical information can be extracted, we analyze color polarization of 
the quark-gluon plasma in a plane-wave Ansatz. 

We start by reviewing the work relevant to hard thermal loops in QCD. The motivation 
that led to their discovery was that physical quantities (such as damping rates) in hot QCD 
were gauge- dependent when computed using the usual loop expansion [0]. The solution to 
this puzzle was first proposed by Pisarski 0. Subsequent development was carried out by 
Braaten and Pisarski 0], and by Frenkel and Taylor These authors realized that, in the 
diagrammatic approach to high-temperature QCD, a resummation procedure is necessary 
in order to take into account consistently all contributions at leading-order in the coupling 
constant. Such contributions were found to arise only from one- loop diagrams with "soft" 
external and "hard" internal momenta. "Soft" denotes a scale ~ gT and "hard" refers to 
one ~ T, where g <^ 1 is the coupling constant, and T denotes the plasma temperature. 
Such diagrams were called "hard thermal loops" in . The HTL approach was successful 
in providing gauge-invariant results for physical quantities. Identifying the momentum scales 
that are relevant to the study of a hot quark-gluon plasma, as was done in P,^, was an 
essential step for all further developments on HTLs. 

An effective action for HTLs was given by Taylor and Wong who, after imposing 
gauge invariance, solved the resulting condition on the generating functional. Efraty and 
Nair have identified this gauge invariance condition with the equation of motion for 
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the topological Cliern- Simons theory at zero temperature, thereby providing a non-thermal 
framework for studying hard thermal physics. Along the same line of research, the eikonal 
for a Chern-Simons theory has been used by Jackiw and Nair P) to obtain a non-Abelian 
generalization of the Kubo formula, which governs, through the current induced by HTLs, 
the response of a hot quark-gluon plasma. 

Another description of hard thermal loops in QCD has been proposed by Blaizot and 
lancu 1^. It is based on a truncation of the Schwinger-Dyson hierarchy and yields quantum 
kinetic equations for the QCD induced color current. These kinetic equations, as well as 
the generating functional for HTLs, were obtained in |0 by performing a consistent expan- 
sion in the coupling constant, which amounts to taking into account the coupling constant 
dependence carried by the space-time derivatives. This dependence is extracted by going 
to a coordinate system which separates long-wavelength, collective excitations carrying soft 
momenta from the typical hard energies of plasma particles. 

Alternatively, Jackiw, Liu and Lucchesi have shown how HTLs can be derived from 
the Cornwall- Jackiw- Tomboulis composite effective action by requiring its stationarity, 
and by using the approximation scheme developed in p. 

The resummation prescription of Braaten-Pisarski and Frenkel- Taylor, as well as the con- 
sistent expansion in the coupling constant developed by Blaizot-Iancu, although remarkably 
insightful, are technically very involved and necessitate lengthy computations. Furthermore, 
they are puzzling with respect to the very nature of hard thermal loops. One wonders if a 
quantum field theoretical description of hard thermal loops (involving gauge-fixing and ghost 
fields) is required. Indeed, we are faced with the following situation: in the resummation 
approach, HTLs emerge from loop diagrams, and in Blaizot and lancu's work, they arise 
from the Schwinger-Dyson equations. However, hard thermal effects are UV-finite since they 
are due exclusively to thermal fluctuations. One might therefore be able to describe such 
effects within a classical, more transparent, context. 

This motivated us to develop a classical formalism for hard thermal loops in QCD, the 
natural starting point being the classical transport theory of plasmas (see for instance [|T^ ). 
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Our effort was encouraged by the fact that for an Abehan plasma of electrons and ions, the 
dielectric tensor computed |T^ from classical transport theory is the same as that extracted 



from the hard thermal corrections to the vacuum polarization tensor Moreover, the 

same situation is encountered for non-Abelian plasmas |1^,|T5|. 

The classical transport theory for non-Abelian plasmas has been established by Elze 



and Heinz [Q, before hard thermal effects were an issue. The HTLs of QCD were not 
uncovered in these early works, mainly due to the lack of a motivation to do so and because 
the transport equations had been linearized, thereby neglecting non-Abelian contributions. 
There has not been, to the best of our knowledge, any attempt to derive the complete set 
of HTLs for QCD from classical transport theory. The aim of the present paper is to give a 
detailed account of this derivation, the results of which have already appeared recently in a 
Letter ||]. 

In our approach, the generating functional of HTLs (with an arbitrary number of soft 
external bosonic legs) arises as a leading-order effect in the coupling constant. We start 



by reviewing the Wong equations for classical colored particles. Following |T^, these are 



substituted into the transport equation, which governs the time evolution of the one-particle 
distribution function, thereby yielding the so-called Boltzmann equation. The latter, aug- 
mented with the Yang-Mills equation relating the field strength to the color current, form 
a consistent set of coupled, gauge-invariant differential equations known as non-Abelian 
Vlasov equations. 

Expanding the distribution function in powers of the coupling constant and considering 
the lowest-order effects, we obtain a constraint on the color current. The latter constraint is 



equivalent to the condition found in |TT|], and previously in 0, on the induced current. The 
constraint on the color current leads to the generating functional of hard thermal loops. 
This work is stuctured as follows. Section || describes classical transport theory for 



the quark-gluon plasma. The latter is reviewed in Subsection |11 A| . In Subsection |IIB| , we 



discuss and justify the ad hoc phase-space integration measure, using Darboux variables. 
Subsection |II Q presents an analysis of the gauge invariance of the system of non-Abelian 



Vlasov equations. Section |TT| contains the derivation of the hard thermal loops of QCD. As 
a consequence of constraints satisfied by the induced current (which are derived in Subsec- 
tion [111 A| ), we obtain the generating functional of hard thermal loops (Subsection [111B| ). In 
Section 0, we compute the polarization tensor from classical transport theory (at leading- 
order in g) and extract the expression for Landau damping. The consistency of our result 
with previous ones is discussed. Section states our conclusions. In particular, we dis- 
cuss there the validity of our approximations. In Appendix 0, we check the validity of the 
Boltzmann equation. Appendix |B| presents a proof of the covariant conservation of the color 
current. 

II. CLASSICAL TRANSPORT THEORY FOR A NON-ABELIAN PLASMA 

A. Classical motion and non-Abelian Vlasov equations 

The classical transport theory for the QCD plasma was developed in [|14|, which we 
follow in this subsection. Consider a particle bearing a non-Abelian SU{N) color charge 
Q", a = 1, ...,N'^ — 1, traversing a worldline x°(r), where r denotes the proper time. The 
dynamical effects of the spin of the particles shall be ignored, as they are typically small. 
The Wong equations \W\ describe the dynamical evolution of the variable^ x^, and Q"". 



m—=p^', (2.1a) 
dr 

^^ = 9Q''FrP. , (2.1b) 

m^ = ~gr'^p^^AlQ^ . (2.1c) 



^Note that we are using the kinetic momentum, rather than the canonical one. A formulation in 
terms of canonical variables would be equivalent IT! 
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The /"^"^ are the structure constants of the group, F^" denotes the field strength, g is the 
couphng constant, and we set c = = fc^ = 1 henceforth. Equation ( |2.1b[ ) is the non- 



Abehan generahzation of the Lorentz force law, and ( |2.1c| ) describes the precession in color 
space of the charge in an external color field AJ^. It is noteworthy that the color charge is 
itself subject to dynamical evolution, a feature which distinguishes the non-Abelian theory 
from electromagnetism. 

The usual phase-space is now enlarged to {x,p,Q) by including into it color de- 

grees of freedom for colored particles. Physical constraints are enforced by inserting delta- 
functions in the phase-space volume element dxdPdQ. The momentum measure 

dP = ^^2e{po)5{p'-m') (2.2) 

guarantees positivity of the energy and on-shell evolution. The color charge measure enforces 
the conservation of the group invariants, e.g., for ^[/(S), 

dQ = d'Q 6{QaQ^ - gs) 6{dabcQ''Q'Q' - gs) , (2.3) 

where the constants q2 and ^3 fix the values of the Casimirs and date are the totally symmetric 
group constants. The color charges which now span the phase-space are dependent variables. 
These can be formally related to a set of independent phase-space Darboux variables. This 
derivation is presented in Subsection |II B| below. 



The one-particle distribution function f{x,p,Q) denotes the probability for finding the 
particle in the state {x,p, Q). It evolves in time via a transport equation, 

m^^^^^^ = C[/](a:,p,g), (2.4) 

where C[/](x,p, Q) denotes the collision integral, which we henceforth set to zero. Using the 
equations of motion ( |2.1| ), ( p.4| ) becomes, in the coUisionless case, the Boltzmann equation: 



00 r\ 

-9QaFl^-gfa,cAlQ' 



f{x,p,Q) = 0. (2.5) 



_^x^^ ''"dp, '^^ dQa_ 

In Appendix ^, an explicit microscopic distribution function is presented and used to check 
the validity of (^). 



A complete, self-consistent set of non-Abelian Vlasov equations for the distribution func- 
tion and the mean color field is obtained by augmenting the Boltzmann equation with the 
Yang-Mills equations: 

= Jf^^ix) . (2.6) 



X] 



The covariant derivative is defined as D"^^ = d^5"'^ + g f A^. The total color current J'^"( 
is given by the sum of all contributions from particle species and helicities, 

J'\x)= E (2-7) 

species helicities 

Each j^"'{x) (species and spin indices are implicit) is computed from the corresponding 
distribution function as 

j>"'{x)=g J dPdQp''Q''fix,p,Q) (2.8) 
and it is covariantly conserved, 

{D,f r (x) = , (2.9) 

as can be checked by using the Boltzmann equation (a detailed proof is presented in Ap- 
pendix For later convenience, we define the total and individual current momentum- 
densities: 

J^^{x,p)=Yl E fi^.P), f\x,p)=g j dQp^Q'^fix.p.Q) . (2.10) 

species helicities 

Note that a solution to the set of Vlasov equations (^]5|)-(^]6|) is specified by giving the 



forms for the gauge potential A^{x) and for the distribution function f{x,p,Q). 



B. Phase-space for colored particles 

In order to carry out the transport theory analysis for classical colored particles, it has 
been necessary to extend phase-space by the addition of the color charges. In (|2.3D, the 
charges are constrained to remain within the group manifold by means of delta-functions 
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which fix the values of the (representation-dependent) group Casimirs. In fact, at an opera- 
tional level, this is the approach adopted in the rest of the present paper. In this subsection, 
we formally justify this approach by analysis of the symplectic structure of the group man- 
ifold We work out exphcitly the SU{2) and SU{3) cases. 

The group SU{2), is generated by three charges, {Qi,Q2,Q3), and has one Casimir, 
Q'^Qa- The structure constants are fahc = ^abc, while dabc = 0. From the point of view 
adopted throughout the rest of this paper the phase-space color measure is 



dQ = dQi dQ2 dQs SiQ'^Qa - q2) 



where q2 denotes the value of the quadratic Casimir. 



(2.11) 



New coordinates (0, vr, J) may be introduced by the following transformation [|T8 



Ql = COS0 V — Tf^ 



Q2 = sin V - vr^ , = n . 



(2.12) 



Note that vr is bounded, — J < vr < J. That the group manifold has spherical geometry 
is readily apparent if one chooses vr = J cos 6. The variables and vr form a canonically 
conjugate pair; the Poisson bracket may be formed in the conventional manner: 

dA dB dA dB 



PB 



(2.13) 



J(f) dir dn d(j) 

It is easily verified that the charges as given by ( p. 12 ) form a representation of SU{2) under 
the Poisson bracket, i.e., 



{Qa, Qb}pB — ^abc Qc ■ 



(2.14) 



The above Poisson bracket structure allows one to identify and vr as Darboux variables 
(see for instance IT^). The Jacobian of the transformation from {Qi,Q2,Q3) to (0, tt, J) 
takes the value 

d{Qi,Q2,Q3) 



J . 



(2.15) 



9(0, TT, J) 

Performing the change of variables ( |2.12D in ( ^.llj ) and substituting the value of the quadratic 
Casimir, Q'^Qa = J^, the color measure reads 



dQ 



'd7idJJS{j'^-q2) 



(2.16) 



which, upon integration over the constrained variable J, is just the proper, canonical volume 
element d(j)d7T, up to an irrelevant constant. 

The group SU{3) has eight charges, {Qi, . . . jQs) and two conserved quantities, the 
quadratic and the cubic Casimirs, Q'^Qa and dabcQ°'Q^Q'^, respectively. The phase-space 
color measure is quoted above in (|2.3|) . 

As in the SU (2) case, new coordinates (0i, 02, 03, tti, 112, tts, Ji, J2) may be introduced by 



means of the following transformations |T^: 



Qi = cos 01 7r_|_ 7r_ , 

Q4 = C++ 7r+ A + C+_ 7r_ 
= C_+ 7r_ A - 7r+ 

in which we have used the definitions: 



Q2 



sm (pi 7r+ 7r_ , 



= 5'++ 7r+ A + S'+_ 7r_ 5 
= S_+7r_ A-S__7r+5 



(2.17) 



7r+ = a/tTs + TTi , 



7r_ 



±± = cos 



i(±0i + V302±- 



5. 



±± 



- TTl , 

sin I (±01 + a/302 ± 03) 



(2.18) 



and A, B are given by 



1 



Ji — J2 



7r3 



7r2 \ / Ji + 2J2 

v^; I 3 



7r3 



7r2 \ / 2 Ji + J2 

v^; I 3 



7r3 



V3y 



1 



27r3\ 



'^2-^1 , vr2\ /J1 + 2J2 



Note that in this representation, the set {Qi,Q2,Q3) forms an SU{2) subgroup with 
quadratic Casimir Q\ + Qi + Q\ = vr|. It can be verified that the expressions above for 
Qi, . . . , Qs form a representation of the group SU{?>): 



{Qai QbjpB — fabcQc 



(2.19) 



under the Poisson bracket 



{Ai?}pB = Ef 

2 = 1 \ 



dAdB OAdB' 



(2.20) 



drii d-Ki d(j)i 
where the canonical pairs are vrj}j=i 2,3. 

As is imphcit in the above, the two Casimirs depend only on Ji and J2. They can be 
computed, using the values given in the table below, as: 

1 



dabcQ'^Q'^Q" = -^{Jl - J2){Jl + 2J2)(2Jl + J2) 
io 



(2.21a) 
(2.21b) 



dabc 


dii% diAQ cii57 d22% 


C?247 (^256 C^338 0^344 C^355 '^366 C^377 


(^448 




d%Q% djYs 


dsss 


Value 


1111 

2 2 V3 


11111 11 
2 2 2 2 2 2 


1 

2v^ 


1 

2^3 


1 1 

2V3 2V3 


1 

V3 



Table: Values of the (non-zero) SU{3) totally symmetric constants. 



The phase-space color measure for SU{3), given in ( |2.3| ), may be transformed to the new 
coordinates through use of ( |2.21| ) and evaluation of the Jacobian 



d{QuQ2,...,Q8) 



9(01, 02, 03, TTi, TTa, TT^, Ji, J2) 

The measure reads: 



V3 
48 



Ji J2 {Ji + J2 



(2.22) 



dQ 



1 

d-Ki dTX2 dn3 dJi dJ2 — Ji J2 (Ji + -h) 5 (-{Ji + J1J2 + Jl) - ^2) x 



5(^(Ji- J2)(Ji + 2J2)(2Ji + J2) 



53 



(2.23) 



Since the two Casimirs are linearly independent, the delta-functions uniquely fix both Ji 
and J2 to be representation-dependent constants. Upon integrating over Ji and J2, (p.23|) 
reduces to a constant times the proper canonical volume element ni=i d(j)id7ii. 

The construction of the canonical phase-space measure for the general SU{N) case is a 
departure from our purposes and will not be undertaken here. Nevertheless, based on the 



examples we treated explicitly, it is apparent that no difficulties will arise for > 3 ||T9|. 

In principle, the classical transport theory analysis should be carried out using canonical, 
independent integration variables and the phase-space volume element should be taken to 



be the proper canonical volume element. In this subsection, we have shown the equivalence 
of the ad hoc phase-space color measure and the proper canonical volume element. Hence, 
the use of the color charges as phase-space coordinates is justified. 



C. Gauge invariance of the non-Abelian Vlasov equations 

Before addressing the question of the gauge invariance of the system of Vlasov equa- 
tions, we consider the Wong equations These are invariant under the finite gauge 
transformation^: 

= x^' , (2.24a) 

= p^ , (2.24b) 

Q = UQU-^ , (2.24c) 

A^ = U A^U-^ --U ^U-\ (2.24d) 

where U{x) = exp[—g e"'{x) t"] is a group element. 

Accordingly, the derivatives appearing in the Boltzmann equation ( |2.5| ) transform as: 

d d 



dxt^ dxi^ 
d d 



2 Iy(| Q]^) (2,25a) 

(2.25b) 



Consequently, the Boltzmann equation (rewritten here in terms of traces): 

f(x,p,Q) = (2.26) 



^We use here matrix notation, e.g. Q = Qat^, qq = qq^ t"", where the generators are represented 
by antihermitian matrices t"' in the fundamental representation, = f"'^'^t^, and they are 

normahzed as Tr (t"t^) = — ^S""^. 
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becomes, in the new coordinates: 



d 



d 



d 



^9P' Tr ([ {g-^U)U-' , Q ]j^)+2 9P^ Tr {QF,,] 



dQ 
d 



_d_ 

dp,. 



/(i,p,Q) =0 



where we have defined 



(2.27) 



f{x,p,Q) = f(^x{x,p,Q),p{x,p,Q),Q{x,p,Q)^ . 



(2.28) 



Simphfying (|2.271) , we obtain: 



f(x,p,Q)=a. (2,29) 



This proves that the Bohzmann equation is invariant under gauge transformations. On 
the other hand, the Yang-Mills equation (|2.6|) is gauge- covariant. Indeed, the color 
current ( p^.8| ) transforms under ( [^.24| ) as a gauge covariant vector: j^{x) — > j^{x) = 
J dPdQp'^ Q f{x,p,Q). Due to the gauge-invariance of the phase-space measure, to the 
transformation property of / (|2.28|), and to (^.24|), j^{x) may be rewritten as: 



fix) 



J dPdQp^'UQU-^ f{x,p,Q) = U f{x) f/"^ . 



(2.30) 



Hence, the system of non-Abelian Vlasov equations is gauge-co variant, with the distribution 
function f{x,p,Q) transforming as a scalar. Note that the gauge symmetry also implies 
that the gauge transform f{x,p, Q)} of a set of solutions {A^{x), f{x,p, Q)} to the 

Vlasov equations: 



A,ix) = UA,ix)U-'-^-U^^U-' 
f{x,p,Q) = fix,p,UQU~') 



(2.31a) 
(2.31b) 



is still a solution. 



III. EMERGENCE OF HARD THERMAL LOOPS 
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A. Constraint on the color current 



Classical transport theory is now employed to study soft excitations in a hot, color- 
neutral quark-gluon plasma. In the high-temperature limit, the masses of the particles can 
be neglected and shall henceforth be assumed to vanish. The wavelength of a soft excitation 
is of order and the coupling constant g is assumed to be small. We then expand the 
distribution function f{x,p,Q) in powers of g: 

/ = /°) + <?/^) + <?¥^) + ... , (3.1) 

where f^^^ is the equilibrium distribution function in the absence of a net color field, and is 
given by: 

f^'\po) = CnBAPo) ■ (3.2) 

Here C is a normalization constant and nB^piPo) = l/(e^'^°' =F 1) is the bosonic, resp. 
fermionic, probability distribution. 

At leading-order in g, the color current ( |2.1CI| ) is 

f%x,p)=g'l dQp>^Q^f'\x,p,Q) , (3.3) 

while the Boltzmann equation ( p.5|) reduces to 



P' ~ 3 /^^nx,P,Q) =p'^g.F;,^/W(po) . (3.4) 



Due to the softness of the excitation, the in the above equation is of order g\A\, so 
we are taking into account consistently all contributions of order g. The approximation 
we use guarantees that the non-Abelian gauge symmetry of the exact Boltzmann equation 
( p.5[ ) is preserved in the approximate equation ( p.4| ). As a consequence f^^^ and f^^\ like 
/, transform separately as gauge- invariant scalars. Other approximations, which have been 



carried out in the past |T^ , have discarded the non-Abelian contributions, thereby breaking 
the non-Abelian gauge symmetry of the Boltzmann equation. 

The equations ( p. 3D and (|3.4| ) yield the following constraint on the color current: 
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[p ■ D nx,p)r = g^v'fFt,^^ i^j dQ Q"Q,/(°)(po)) , (3.5) 

where, from color symmetry, we have/ dQ Q°'Qbf''^\po) = Cb,f ^b,f{Po) ^1 with Cb = 
N, Ci? = I for gluons, resp. fermions. Thus, upon summation over all species {Np quarks, 
Np antiquarks and one [(A^^ — l)-plet] gluon) and helicities (2 for quarks-antiquarks and for 
the massless gluon), (|3.5|) yields. 



[p-D J^^ix,p)r = 2g^p'^p''F:,-^[NnB{po) + NFnF{po)] . (3.6) 

dpo 



Similar results have been obtained in PJTl[], in a quantum field theoretic setting. 



B. Derivation of hard thermal loops 

Subsequent steps which lead to the generating functional of HTLs have been described 



m 



IT]]] , the results of which were used straightforwardly in [|[], for the sake of brevity. Here, 
we present a simpler derivation of HTLs by exploiting fully the structure of the momentum 
integration measure ( p.2|) . 



We first integrate equation ( p.6| ) over |p| and Pq using the massless limit of the momentum 



measure dP (|2.2|) . Therefore, the (massless) mass-shell constraint enforces |p| = Po, and we 
thus introduce the unit vector p = p/|p|. Introducing also v = (l,p), the integration of 
( p.6| ) yields (group indices are henceforth omitted): 



vD J^{x,v) = -2TT^mlv''vf Fpo{x) , (3.7) 
where rriD is the Debye screening mass 



N + Nf/2 

'mD = 9Td , (3.8) 

and we have defined 

J''{x,v) = J \p\^d\p\dpo2e{po)6{p^) J^{x,p) . (3.9) 
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Notice (for later use) that, using / dP = J j^^\p\'^d\p\dpo29{po) 5{p^), where d^l denotes 
integration over all angular directions of the unit vector p, we can rewrite the expression 
J^{x) = J dPJ^{x,p) for the color current as 

dVL 



After decomposing J^{x,v) as 



(27r)= 



(3.10) 



(3.11) 



we get as our final condition on the color current: 



d 



V ■ D J^(x, v) =2 71^ ml v^"— (v ■ A(x)) 

ox" V / 



(3.12) 



It has been shown that solutions to (|3.12|) can be obtained from a functional W^(v4, v) 
as H 

SW{A,v) 



6A^{x) 



(3.13) 



Equation ( p.l2|) then implies that W{A,v) depends only on Aj^ = v ■ A, i.e. W{A,v) 
W{A+), and J'^ = ^^^'^^ v^. In turn, W{A+) satisfies, as a consequence of ( p. 121) , 



(3.14) 



By introducing new coordinates ( 



X^ = V ■ X , X_ = V ■ X 



x± = X - (p ■ x)p , 



(3.15) 



with V = (1, — p) and x_l ■ p = 0, we can rewrite f ■ ^ as 9+ and ( |3.14| ) becomes: 



A, 



6W{A^ 
~6AZ 



2r:^ml±,A,. 



(3.16) 



Now using ( p.lO| ), ( |3.11| ) and ( |3.13| ), we define an effective action T that generates the 
color current, i.e., J^ix) - 



' Xi^lx)^ ' where V takes the form: 



r = !^ I d^xAl{x)Al{x)- ^ 



W{A^ 



(3.17) 
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This is the expression for the effective action generating hard thermal loops while 
equation ( |3.16| ) represents the condition of gauge invariance for this generating functional. 
By solving (|3.16| ) , Taylor and Wong , as well as Efraty and Nair [j^ , have given an explicit 
form for the functional in the second term of ( p.l7| ). The first term is a mass term 

for Aq{x) and describes Debye screening. 

This concludes our derivation of the hard thermal loops of QCD from classical transport 
theory. 

IV. APPLICATION: COLOR POLARIZATION 

As an application of the classical transport formalism presented above, we solve the 



approximate Boltzmann equation ( p.4|) for plane-wave excitations in a collisionless isotropic 
plasma of quarks and gluons. Recall that in the case of a collisionless plasma of electrons 
and ions the Abelian version of equation ( |3.4| ) has been solved exactly for an electromagnetic 
plane- wave , making it possible to study the response of an Abelian plasma to a weak 

field. We shall proceed analogously in the non-Abelian case. We consider a plane-wave 
Ansatz in which the vector gauge fields only depend on through the combination x ■ k, 
where k'^ = (c(j,k) is the wave vector, i.e., ^^(x) = ■ x). With this Ansatz (which has 
been used in |jT^ to study the non-Abelian Kubo equation) the solution of (|3.4|) is 



/«(a;,p,g) =g„ (A^,{x)-u;^^^j ^/'Hpo) • (4.1) 
Hence, the color current is given by 

The integration over color charges can be done by using / dQ Q'^Qb = Cb^f^I with Cb = 
N, Cp = \ for gluons, resp. fermions. The integration over po and |p| is straighforward 
as well. Upon summation over all species and helicities (see Section III for notations and 
conventions), we get the following expression for the total color current: 
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V^UJ 



V ■ AJx) 



Att \ V ■ k 

The polarization tensor li^ can be computed from ([4 .31) by using the relation 



(4.3) 



It reads: 



= j d'yIiZ{x - y) Al{y) . 



Klix ~y) = ml {-g^'g'^' + /'^^(o;, k)) ^^(x - y) S^t 



(4.4) 



(4.5) 



where I'^'^ is defined as 



(4.6) 



— k ■ V 

To avoid the poles in the above integrand, we impose retarded boundary conditions, i.e., we 
replace u hj u + ie. Using the identity 



1 1 

= V in S(z) 

z + ie z ^ ^ 



(4.7) 



where V stands for the principal value, the real and imaginary parts of the polarization 
tensor are 



Ren::,^(^,k) = -6., ml (-^^V^P / S ZT^^T^^ 
Imn(^^((x;,k) = -5ab mlnu J ^ v'' S{uj - k ■ v) . 



(4.8a) 
(4.8b) 



The imaginary part of the polarization tensor (|4.8b|) describes Landau damping in the 
quark-gluon plasma. Explicitly: 



Imn°'^(u;,k) = -6abmln 



Im n^^b(^> = -6r,h mn vr 



CO 



2|k| 

2 /'ik 



UJ 



A\k]^[~']k\ \ ~Wr2\k\^ ik|2 



(4.9a) 
(4.9b) 

\k\^-uj^) . (4.9c) 



From the above 6'-functions it is apparent that Landau damping only occurs for color fields 
with space-like wave vectors. This is also true for an Abelian plasma 
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Evaluation of the real part of the polarization tensor ( |4.8a|) yields 



Ren^,(cu,k) = 5, 



ab 



where 



Yii{u,k) = m\) 



ijj 
%k\ 



In 





\k\ 


Ijj — 


\k\ 



nf(u;,k) = -m^ 



2|k| 



1 + 



|k| 

UJ 



UJ 



In 





\k\ 


UJ — 


\k\ 



(4.10a) 
(4.10b) 
(4.10c) 

(4.11a) 
(4.11b) 



Equations ( [4.10 )-( 4.11 ) characterize Debye screening, as well as longitudinal and transverse 
plasma waves. 

Our results for the HTLs of the polarization tensor agree with those obtained in the high 
temperature limit using quantum field theoretic techniques^ ||T5|J^,p|,^p|,|r?| . We emphasize 



that the above results are gauge-independent, and obey the Ward identity 



(4.12) 



as should be expected from the gauge invariance of our formalism. 

Previous applications of classical transport theory to QCD have utilized an Abelian- 
dominance approximation to compute the polarization tensor ||14|. It is noteworthy that 
there one recovers the same values of the polarization tensor that we found here. The 
reason for this agreement is that the leading-order contribution to the color current is made 
linear in the gauge field by the plane- wave Ansatz [|1^], exactly as happens in the Abelian- 
dominance approximation. However, the Abelian-dominance approximation cannot give a 
proper account of the whole set of HTLs, such as thermal corrections to n-point functions, 
n>3. 



^The connection between the retarded polarization tensor computed here and the time-ordered 
polarization tensor that is commonly used in quantum field theory has been studied in |^. 
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V. CONCLUSIONS 



In this paper, we have shown how classical transport theory can be used to derive the 
hard thermal loops of QCD. This formalism, we believe, is more direct and transparent than 
previous approaches based on perturbative quantum field theory. Indeed, hard thermal 
loops represent UV-finite thermal corrections to propagators and vertices. They arise from 
thermal scattering within a hot assembly of particles, and one may reasonably expect them 
to be describable in terms of classical physics. 

The fact that we are modeUing the high-temperature, deconfined, phase of QCD allows us 
to treat color classically, and enables the colored constituents of the plasma to be identified as 
quarks and gluons. Employing classical transport theory to study colored particles requires 
incorporating the color degrees of freedom into phase-space. A consistent measure must be 
defined over the new color coordinates. Furthermore, conservation of the group Casimirs 
under the dynamical evolution must be ensured. One means to accomplish these goals is to 
include delta-function constraints into the phase-space volume element. We have formally 
justified this ad hoc procedure by relating the dependent color degrees of freedom to a set 
of independent Darboux canonical variables, and by proving that the corresponding volume 
elements are equivalent. 

A system of non-Abelian Vlasov equations describes transport phenomena in the QCD 
plasma. This system governs the evolution of both the single-particle phase-space distribu- 
tion functions and the mean color fields. It would be a formidable task to solve the transport 
equation in the most general case, hence suitable approximations must be made. 

First, we specialize to a coUisionless plasma, in which there is no direct scattering between 
particles. This situation is not devoid of interest since collective mean-field effects can, and 
indeed do, arise. 

Second, we employ a perturbative approximation scheme. We assume that the plasma 
is near equilibrium and expand the phase-space distribution function in powers of the 
gauge coupling constant. At the high temperature which must prevail for the formation 
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of a quark-gluon plasma, the thermal energies of the particles are sufficiently large that 
the effects of the interactions with the gauge fields are comparatively small, and we expect 
perturbation theory to be valid. 

Taking the high temperature limit constitutes our third approximation. In this limit, the 
masses of the particles can be neglected. Furthermore, the plasma is in a highly degenerate 
state, so that the equilibrium distribution functions are determined by the spin-statistics 
theorem. 

We demand that gauge invariance be preserved by our perturbative expansion. In its 
lowest, non-trivial, order this expansion leads to the generating functional of HTLs. That 
gauge invariance is the appropriate guiding principle in uncovering hard thermal loops is 
not surprising. To apply this principle to the quantum field theoretic calculation of HTLs 
involves gauge fixing, ghosts, and resummation of classes of Feynman diagrams. In contrast, 
the route that one has to follow in order to adhere to the gauge principle is straightforward 
within classical physics. Therefore, the gauge invariance property of hard thermal loops is 
self-evident in our formalism. 

Acknowledgements: We thank Professor R. Jackiw for constant encouragement and 
for many useful and enjoyable discussions. 

APPENDIX A: MICROSCOPIC DESCRIPTION 

In a microscopic description, the particle's trajectory in phase-space is known exactly. 
With this knowledge, we can construct a distribution function f{x,p,Q) (without loss of 
generality, we shall consider only one particle): 

f{x,p,Q) = [- S^'\x-x{t)) 5W(p-p(r)) 5(^'-^)(g-g(r)) , (Al) 
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where x{t), p{t) and Q{t) obey the Wong equations (|2.1| ), i.e. they naturally fulfill the 
mass-shell and Casimir constraints. For convenience, those constraints are here subsumed 
into the distribution function instead of being contained in the phase-space volume element. 
Had we not done this, ( |Al| ) would have to be written in terms of both a 3-dimensional 
5-function in momentum space, and an A^(A^ — l)-dimensional 5-function in color space. 
We now prove that the expression (|A1|) for / satisfies the collisionless Boltzmann equation 



]5|). The first term in ( p.5|) can be rewritten, by using the properties of the (5-function, as: 



d 



/ = - 



dr 



m 



d 



X — X(T 



5W(p-p(r)) 5(^^-1) (g-Q(r)) . (A2) 



^x^'{T) 

From this, after using the Wong equation for the variable x^{t) and applying the chain rule, 
we get: 



P'£-J=-Jdr ^ {x - x{t)) 5W (p - p(r)) 6^^'~'^ (Q - Q(r)) 
Similar arguments yield, for the second term in the Boltzmann equation: 

- 9P'QaF;;^^f = -J rfr5W(x - x{t)) \^S^'^{p - p(r))] d^'^'-'^Q - Q{ 



(A3) 



(A4) 



and for the third term: 
d 



(A5) 

Adding together the equations (|A3|) , (|A4|) and ([A5|), we observe that the left hand side - 
the r-integral of a total r-derivative - vanishes, thereby yielding the collisionless Boltzmann 
equation ( |2.5| ). 



APPENDIX B: CONSERVATION OF THE COLOR CURRENT 



Let us verify that the color current ( |2.8| ) is covariantly conserved. Using the collisionless 
transport equation, one can compute 
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dt,f''{x)^gj dPdQp''Q''dJ{x,p,Q) 



^g'l dPdQp^Q^ |g,F^^^(^)^ + /'^fc^^,(^)g^_^j f{x,p,Q) , (Bl) 



d 



where the color measure is dQ — d^^^ ^^Q C{Q), and C{Q) specifies the color constraints 
in phase-space. Integrating by parts and discarding surface terms, one gets 

d^f'ix) = JdP I dQ{Q'Q,gi^F'^,ix) + p^5,df"'A,,ix)Q' + p^f''^A^,{x)52Q:) 



+ 



I S^'-'^Q p^A,,{x)f"^Q,-^C{Q) 



f{x,p,Q) . 



(B2) 



Among the four terms in the right side, only the third one survives. The first two terms 
cancel due to antisymmetry of F^^ and f^"^, respectively. The last term also cancels, since 
the constraints C{Q) are gauge-invariant, i.e., 



(B3) 



where 5(5" denotes an infinitesimal gauge transformation with arbitrary parameter eh{x). 
(For SU{3) this last property can be explicitly checked by using the Jacobi-hke identity 
fabcddec + fadcdebc + faecdbdc = 0.) Finally, ouc obtaius the expression for the covariant 
conservation of the color current: d^j^°'{x) + 5'/"'"^^^6(a;)j^(a;) = 0. 
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